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Abstract 

We propose an "interferometric" scheme for Bose-Einstein condensates using 
near-held diffraction. The scheme is based on the phenomenon of intermode 
traces or quantum carpets; we show how it may be used in the detection of 
weak forces. 
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One of the most important aspects of the study of Bose-Einstein Condensates (BEC's) 
in dilute trapped alkali gases has been the development of the matter-wave equivalent of 



accompanied by that of useful atom optical elements to manipulate the flowing condensates. 
We have proposed previously a method to split and focus a condensate using an external 
potential [[|. In addition to direct focusing and beam splitting, an integral part of the field 
of atom optics is matter- wave interferometry ||. In this Brief Report, we propose the use 
of what is essentially single slit diffraction for matter-waves to perform interferometry with 
BEC's which should be realisable with the current experimental technology Several theoret- 
ical studies of such phenomena for ordinary non-condensed wave packets have been carried 
out recently |4j], showing in the spacetime probability distribution the so-called "quantum 
carpets." The shapes can be explained by a careful consideration of the solutions of appro- 
priate Schrodinger's equation. With the uncondensed atoms, however, at most few atomic 
wave packets are allowed at a given time in order to maintain the quantum coherence re- 
quired for the spatio-temporal patterns. Clearly this poses a serious experimental challenge; 
for instance, an exceedingly good vacuum system would be needed to maintain the carpet 
for a reasonable duration of time. An immediate advantage of BEC is that the patterns 
can be formed with N atoms, which allows possible experimental observation of the atomic 
quantum carpets, and its application as a matter-wave interferometer. 

In order to construct a quantum carpet, we begin with a Bose-Einstein condensate close 
to temperature T = 0, stored in a one dimensional harmonic trap of frequency u. (The 
minimum of the trap is set to be at x = 0.) At time t = we suddenly switch off the 
trap and release the condensate in a "harder" potential of the form where £ 3> 2, with 
the walls symmetrically located at x — ±L. For simplicity, we propose to use an inverted 
super-Gaussian, Vb ox ( r ) as our potential: 



where L is defined to be half of the total transverse length of the potential in harmonic 
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optical lasers, the atom laser [![]. It is essential that the development of an atom laser be 
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oscillator units of length, and a = 0.9L. With m = 42, the potential closely approximates a 
box of length 2L. Experimentally, the potential is expected to be well approximated by using 
parallel sheets of optical fields. It is clear that Vb ox is equivalent in the optical context to 
imposing a slit, which simply provides a set of boundary conditions for Maxwell's equations. 
We note that the group of Ertmer has already implemented coherent reflection of BEC from 
an atom optical mirror and also investigated BEC in an atom-optical wave guide, which are 
very closely related to the experimental technology required to observe a quantum carpet 

We describe the dynamical evolution of Bose-Einstein condensates using the Gross- 
Pitaevskii Equation (GPE), which provides an excellent description of BEC at temperatures 
T < T c where T c is the critical temperature of phase transition as confirmed by various ex- 
periments to date. This equation has the form: 
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with C = NUq = 47Th ^ [ Na , where iV is the number of particles, a the interatomic s-wave 
scattering length and M the mass of a single particle. V(r) is the trapping potential which 
may be switched off and replaced by Vb ox (r). The prepared condensate wave packet of 
parabolic shape is evolved according to the one dimensional version of Gross-Pitaevskii 
equation: 
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Here the size, k = C/^o^q) of the nonlinearity is a reduction from the three- to the one- 
dimensional case. Such ID simulation does indeed have a relevance in reality. For example, 
a 3D trap with a box potential in one direction, along with tightly confining harmonic 
potentials in the remaining two directions may be used. Assuming weak directional coupling 
in the condensate motion, this would then produce interference in one direction only, giving 
results effectively for a ID carpet. 

As a result of switching to V^ ox , the kinetic energy term in GPE dominates and the 
condensate wave packet rapidly spreads as function of time, and reflects at the walls acting 



as atomic mirrors. It then interferes with itself with varying degree over time producing 
complex intensity profiles or "fringe patterns." Shortly afterwards when these fringes are 
viewed in the time domain, regular dark canals of low density are formed and we observe a 
remarkable emergence of spatio-temporal patterns 0. The pattern for L = 25aHo ; T = 
and evolution time of to IOtho where ano is one harmonic oscillator length and tho is 
one period of the harmonic oscillator is shown in Fig. [j]. The total duration of simulation 
is therefore approximately 100 ms for a trap with frequency u = 2n x 100Hz. For this 
simulation, we used a nonlinearity k = 500 hu^ano which corresponds to approximately 
20000 Rubidium atoms. Qualitatively similar canals and ridges occur when we propagate 
the same wave packet according to the GPE with k = (i.e. linear, uncondensed case). 
However three fundamental differences occur in the carpet with the nonlinearity turned on: 
(i) The initial spreading is much faster; this can be related to the time dependent excitation 
of higher momenta due to interactions, (ii) The canals and ridges develop a substructure 
along the straight lines predicted by the linear propagation. (Hi) We note complicated 
avoided crossings when two canals cross each other. The most important difference is that 
even for much longer integration times the fractional revivals of the wave function do not 
appear as seen in the linear case. We note that the sizes and the resolution of the observed 
patterns depend strongly on factors such as the strength of nonlinearity k and the length 
2L of the box. In particular, the pattern becomes more distinct and the canals turn to 
ridges and vice versa with fewer number of atoms or smaller k. From an experimental point 
of view, we note that sufficient temporal resolution to observe the carpet structure exists 
with, for instance, the nondestructive imaging procedure which has the shortest probe pulse 
duration of order 1 ms ||. In particular, with a set of repeated measurements made at 
various different reference times one may expect to improve the resolution further. 

We note that the quantum carpets woven with uncondensed atomic wave packets have 
been extensively analysed using various methods. They are found to exhibit regular spatio- 
temporal patterns consisting of straight lines of varying gradients, representing harmonics 
of a fundamental wave velocity Vf = uif/kf = tnr/(2mL), where uuf = hir 2 /(2ML 2 ) and 
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kf = tt/L denote the fundamental frequency and wave number in a box potential of length 
L. The probability density results from the interferences between degenerate eigenmodes 
of the system. The straight lines in the carpet formed from uncondensed wave packets 
represent wave velocities in integer multiples of Vf, vi n = (I ± n)Vf where I and n are 
integers, I, n = 0, 1, . . .. One expects similar analysis for the quantum carpets with BEC, 
and indeed as k — > this is precisely what happens with the BEC quantum carpets. 

However, for large values of k, the self-interaction present in BEC considerably com- 
plicates the analysis of the carpet, leaving numerical analysis as essentially the only viable 
option. However, it is found that one can gain useful insight into the existing structures by 
combining two facts about BEC. One is that, depending on the sign of the s-wave scattering 
length, condensates are able to support bright, grey and dark solitons which are the unique 
features of nonlinear wave equations such as the GPE. In fact, under GPE, any solution 
which is not an eigenstate of the given external potential has a tendency to form solitons 
as t —>■ oo. In particular, a study of grey soliton formation in a condensate resulting from 
collisions have been carried out by Scott et al. |5J Grey solitons, which are the localised 
reduction in the ambient density, or the "density notch" solitons, are found in inhomoge- 
neous condensates with positive s-wave scattering lengths, and they are found to play an 
important role in governing the bulk motion of the condensate. They tend to store energy 
in a slow moving structure and act as a throttle to retard the movement of the conden- 
sate, possibly playing an important role in dissipating superfluid flow ||. Their speed of 
propagation v is controlled by the difference in phase 5 between the parts of the condensate 
separated by the density dip, in a form reminiscent of Josephson effect jTIJ: v = cb sin(5/2) 



where cb = \j Antfap/m 2 is the Bogoliubov speed of sound with p being the condensate 
density. A corollary is that the speed of grey solitons provides a sensitive measure of the 
phase difference between adjoining regions of the condensate. 

The second fact comes from the work on soliton dynamics in BEC. Kosevish [|11|] and also 



Reinhardt and Clark [|ItJ have shown that, in the case where Vt rap varies slowly on the scale 
of the condensate correlation length, the grey solitons behaves as a "classical" particle of 



positive mass. The position of the center of the soliton x c (t) obeys the Newtonian equation 
of motion: 

m = (4, 

This has been demonstrated using V tr:ap of harmonic form, which clearly shows the soliton in 
a simple harmonic motion. The main point to note is that the solitons exhibit particle-like 
behaviour, demonstrating, quite remarkably, the possibility of applying classical Newtonian 
mechanics in the description of a purely quantum system where wave-mechanics is supposed 
to dominate its behaviour. 

Ignoring the more fine lines (smaller solitons) one can see immediately about 30 main dark 
straight lines in Fig. p] criss-crossing each other (V tTSu p = for —L < x < L). They represent 
grey solitons bounced off the side walls, each travelling at some constant speed, determined 
by the amount of phase slip the soliton gains as a result of interference that occurs when the 
condensates reflect from the walls. The different velocities may be understood from the fact 
that upon reflecting and interfering, the condensate gains phase gradient which is step-like, 
yet with different gradient per step as shown in Scott et al. ||. One can see that as soon as 
the condensate meets the wall at time ~ 1tt/u, it creates (or "launches") approximately 15 
major grey solitons travelling at different constant speeds from either side. Faster solitons 
move further away from the walls within the same time duration than the slower ones, and 
reach the opposite wall earlier. The largest velocity is approximately that of the Bogoliubov 
speed in average atomic density. Once reaching the opposite wall the solitons rebound but 
simultaneously generate even more solitons of smaller sizes. It is also noted that the smaller 
grey solitons represented by the finer lines on the carpet continually collide and interact 
with other solitons to eventually alter the paths, sizes and shapes of the major grey solitons 
quite significantly from the simpler Newtonian picture. However, initially at least i.e. up 
to time of ~ 15tt/u, the main grey solitons can clearly be observed to obey the Newtonian 
mechanics. The avoided crossings are indicative of the repulsive interaction. 

We now illustrate the effect of adding a linear variation V to the inverted super-Gaussian 
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potential of the quantum carpet. Such a potential models the effect of an external force, 
and subsequent phase shift across the quantum system. Expected Newtonian trajectories for 
particles travelling at different but evenly distributed initial speeds under several different 
constant accelerations are given in Fig. |2|. Fig. shows the BEC quantum carpet where 
V = 0.05x in scaled units and Fig. ^ the case where V = O.lx in scaled units. In real units 
(in the case the trap frequency is 2n x 100 Hz and using Rubidium atoms) these correspond to 
0.003(7 an d 0.006g respectively where g is the acceleration due to gravity. It is quite clear that 
the canals show exactly what one expects from standard kinematics. This result shows that 
one can deduce the nature and magnitude of the potential (or equivalently interferometric 
phase shifts) from observing the motion of the grey solitons (or the interference fringes). 
This allows not only the ability to detect standard time- independent linear potential but, 
in principle, also gives us the ability to deduce time- varying potential of different shapes, as 
long as the corresponding Newtonian motion can be identified. 

The limitation of any interferometric method would be the intrinsic phase diffusion within 
the condensate which may eventually destroy the interference patterns. In real experiments, 
the phase diffusion rate is expected not to differ too greatly from that of the original MIT 
interference experiment []12||, which reported a decoherence time of the order of 100 ms. 



Assuming that the three-body recombination rate is more or less similar for different atomic 
species, the diffusion time scale is clearly much longer than the integration times we have 
presented here, reflecting the surprising robustness of the BEC interference fringes against 
decoherence. We note also that recently Burger et al. has observed the evolution of dark 
solitons in BEC, created using a phase imprinting method |13|]. Obviously their experiment 



is not a quantum carpet phenomenon; however it is clear that the force measurements 
discussed in this paper may directly be implemented in their experiment. 

In summary, we have studied how one may use the quantum carpets with BEC as an 
interferometer. The spatio-temporal patterns could be explained using the idea of soliton 
formation, along with the classical equation of motion for the solitons. This is a very 
different kind of interferometer than the standard matter wave interferometer such as the 



Mach-Zehnder, and is of interest from a number of different aspects. 
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FIGURES 
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FIG. 1. Quantum carpet with BEC, represented as spacetime probability density. The time and position 
are given in units of harmonic oscillator (h.o.) period and length respectively. 
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FIG. 2. Newtonian trajectory of classical particles under the potentials of the form (a) V{x) = (b) 
V(x) = O.Olx (c) V(x) = 0.05x (d) V(x) = 0.1a;, where x is along the vertical axis. The different lines 
correspond to different initial velocities. 
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FIG. 3. Quantum carpet with BEC, in an externally imposed potential of the form V(x) = 0.05x in 
harmonic oscillator (h.o.) units. 
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FIG. 4. Quantum carpet with BEC, in an external potential of the form V(x) — O.lx. 
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